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Abstract 

An atomistic approach is used to investigate finite-temperature properties of ferroelectric nan- 
odots that are embedded in a polarizable medium. Different phases are predicted, depending on the 
ferroelectric strengths of the material constituting the dot and of the system forming the medium. 
In particular, novel states, exhibiting a coexistence between two kinds of order parameters or pos- 
sessing a peculiar order between dipole vortices of adjacent dots, are discovered. We also discuss 
the origins of these phases, e.g. depolarizing fields and medium-driven interactions between dots. 
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Ferroelectric nanostructures (FENs) are of current high interest, because of their tech- 
nological promise in leading toward miniaturized devices^, and because of their potential in 
resulting in new phenomena. For instance, is was recently predicted that isolated nanodots 
of ferroelectrics can have a vortex structure for their dipoles below a critical temperature 2 ^*^. 
Such vortex does not create any polarization but rather generates a macroscopic toroidal 
moment, that involves the cross product between the vectors locating the i unit cells and 
their local electrical dipoles pj - i.e., it is defined as G = ^ J2i r i x Pi; where the sum runs 
over the N unit cells of the nanodot. Being able to switch the direction of G opens exciting 
opportunities for nanomemory devices 2 ' 5 . 

Interestingly, the possibility for the dipoles to form a vortex structure has been omitted 
in the previous studies (see, e.g., Refs.^»^ that used an effective medium approximation) 
aimed in determining the properties of ferroelectric or paraelectric nanoparticles that are 
embedded in a polarizable medium. As a result, original features may have been overlooked 
in these composite systems. For instance, one may wonder if novel states occur in such two- 
component materials. In particular, can dipole vortices in dots immersed in a polarizable 
medium coexist with a spontaneous polarization (unlike in isolated dots^)? Similarly, 
how do toroidal moments of neighboring dots organize themselves? (e.g, are they lying 
along parallel or antiparallel directions?). It would also be worthwhile to determine the 
dependency of these novel states (if any) on the ferroelectrics strengths of the materials 
forming the dot and medium, and to reveal their governing mechanisms. The aim of this 
Letter is to address the issues mentioned above, by performing atomistic simulations. 

Our atomistic scheme is based on the construction of an effective Hamiltonian, with the 
total energy E being written as a sum of two main terms, 

Ed^iiVHjAvAiWj}) = £avc({u t },{??F}, {??/}) + £loc({u 8 }, {??/}), (1) 

where E &ve - as in Refs.£±a - is the total energy associated with the hypothetical simple A< 
B >C>3 system resulting from the use of the virtual crystal approximation (VCA)^i to mimic 
A(B',B")03 compounds; E\ oc gathers alloying terms going beyond the VCA approximation; 
Uj is the local soft mode (directly related to the electrical dipole) in unit cell i; {t]h} and {rjj} 
are the homogeneous and inhomogeneous strain tensor, respectively 12 ; {o~j} characterizes 
the atomic configuration, that is, = 1 (resp. —1) if there is a B' (resp. B") cation at 
the B-lattice site j of the A(B',B")03 materials. For E ave , we use the analytical expression 



2 



proposed in Ref.— for simple ABO3 systems. For E\ oc , we use the following expression: 

£loc({U;}, {777}, {<Tj}) = ^2[Qj,iCTj Gji ■ Ui + Rjjffj fji ■ Vi] + Are(<7i) u ] ( 2 ) 

ij i 

where the sums over i and j run over unit cells and mixed sublattice sites, respectively, {vj} 
are dimensionless local displacements which are related to the inhomogeneous strain variables 
inside each celi^. Qj,i{<Jj) and Rj ti ((jj) characterize the strengths of the alloying-induced 
intersite interactions, is a unit vector joining the site j to the center of the soft-mode 
vector Uj, and is a unit vector joining the site j to the origin of the displacement v*. 
Practically, we included contributions up to the third neighbors for Qj,i(<Jj), and over the 
first-neighbor shell for Rj^tjj). Note that, for systems made of AB"03 dots embedded in 
a AP/O3 medium, these contributions only play a role near the dots' surfaces - because 
of the analytical expression of the first two terms of Eq. (2). The last term of Eq. (2), 
which involves the Are(erj) parameters, characterizes the onsite contribution of alloying. It 
is an original contribution since it was not included in the alloy effective Hamiltonians of 
Refs.^iH. We incorporate such term here - which is consistent with first-principles results 
of Ref.— - because it provides an easy way to artificially adjust the "ferroelectric strength" 
of AB'03 and AB"03 simple systems, by playing with the An(<Ji = +1) and Aft(<7j = —1) 
parameters. For instance, a large negative (respectively, positive) Are(f7j = —1) leads to a 
strong ferroelectric instability (respectively, no ferroelectric instability) of the pure AB"03 
material. In the following, we will denote Are as the difference between the two alloying- 
onsite parameters, i.e. An = Aft(crj = +1) — Are(o"j = —1). At the exceptions of the 
adjustable Are(<7j = +1) and Are(o"j = —1) variables, all the parameters of our toy model 
described by Eqs. (l)-(2) are those derived for the Pb(Zr .5Ti 5 )0 3 solid solution from first- 
principles calculations 9 . [Such parameters yield a tetragonal ferroelectric ground-state with 
a polarization pointing along a < 001 > direction and a Curie temperature ~ 1000 K, in the 
Pb(Zr .5Ti .5)O 3 bulk]. As a result, we numerically found that Ak(cTj = +1) or Are(<7j = —1) 
~ +0.0094 a.u. is the highest value to have a ferroelectric ground-state in pure AB'03 or 
AB"03 material, respectively. 

The total energy of Eq. (1) is used in Monte-Carlo simulations to obtain the local mode 
vectors, the toroidal moment of polarization and the homogeneous strain tensor {rjjj} at 
different temperatures. Typically, 2xl0 4 Monte-Carlo sweeps are first performed to equili- 
brate the system, and then 8xl0 4 sweeps are used to get various statistical averages. The 
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temperature is decreased in small steps down to 1 K. We also heat up the investigated sys- 
tems by small temperature steps from 1 K, to look for thermal hysteresis. We also computed 
the dielectric susceptibility, as well as the toroidal susceptibility (which characterizes the 
response of the toroidal moment to the curl of the electric fielcU^), as in Ref.— . 

Figure 1 shows the temperature-versus-A/? phase diagram of a 16 x 16 x 16 periodic 
supercell (20480 atoms, 64A lateral size) containing an AB"03 cubic dot of 48A of lateral 
dimension embedded in a host matrix made of pure AB'03. Practically, for positive Ak, 
An(ai = —1) is set to zero while Aftf^ = +1) is allowed to vary. These cases thus correspond 
to a specific ferroelectric dot immersed in different media that are all ferroelectrically-harder 
than the dot and that can either be ferroelectric (small positive Ak) or paraelectric (larger 
positive Ak). The reverse situation applies for the case of negative Ak: Artier, = —1) is 
varied while AK(<7i = +1) vanishes, implying that we mimic either ferroelectric dots with 
smaller ferroelectricity strength than the medium (for small negative Ak) or paraelectric dots 
(for larger negative Ak) inserted in a matrix made from a specific ferroelectric material. 

Figure 1 reveals the existence of six different phases, for which the associated dipole pat- 
terns are displayed in insets [Phase boundaries were practically determined by identifying 
the peak or abrupt jumps of the susceptibilites altogether with the appearance of the spon- 
taneous polarization or toroidal moment. Note that, in principle, further study could reveal 
some of these features to be associated with continuous changes in configuration rather than 
true phase boundaries per se. However, some preliminary results obtained by us on the basis 
of the finite-size scaling analysis 16 do confirm the existence of true phase transitions.] Two 
of these phases are expected based on previous knowledge of ferroelectrics: the paraelectric, 
PE, state at high temperature, and the ferroelectric, FE1, phase occurring at intermedi- 
ate and low temperature when Ak is small in magnitude. In this FE1 state, both the dot 
and medium develop homogeneous parallel dipoles - with these dipoles being larger (resp., 
smaller) in the dot than in the host when Ak is positive (resp., negative), as consistent with 
our definition of Ak (i.e., the difference between the onsite parameters of the medium and 
the dot). 

Four phases of Fig. 1 can be considered as novel structures. They are denoted as FT, 
FE1+FT, FE2 and FE3, respectively. The FT phase, which occurs for the largest positive 
Ak values, exhibits a finite toroidal moment. The appearance of this state in this region 
of the phase diagram results from the fighting of the dipoles in the dot against the large 
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enough depolarizing field^^ (which arises from the significant non-similarity between the 
"more ferroelectric" dot and '"less ferroelectric" host material). The present discovery of 
this FT phase is, in fact, consistent with the previous finding that isolated ferroelectric 
dots surrounded by vacuum exhibit a vortex structure for their dipoles below a critical 
temperatur o 2 i 3 i 4 , because one can think of vacuum as a medium having an infinite positive 
value for An(ai = +1). However and unlike in the vacuum, the dipoles of the host matrix 
in the FT phase become slightly polarized by the nearby dipoles located inside the dots and 
near the surfaces, when this host matrix is still rather "soft" (see corresponding inset of 
Fig.l). For such cases, the medium thus also generates a toroidal moment that is parallel 
to, but of lower magnitude than, the one solely associated with the dot. 

Moreover, the FE1+FT phase appearing in the phase diagram at small temperature and 
within a narrow range of positive Ak is rather remarkable because it displays an unusual 
cohabitation between the toroidal moment and the spontaneous polarization. In this phase, 
the medium generates an electric field inside the dot, that (as in some magnetic nanodots 
under an external magnetic field^ or in some ferroelectric dots subject to electric fields^) 
leads to the shift of the vortex structure with respect to the center of the dot (see the cor- 
responding inset in Fig. 1) and thus activates a polarization. Note, too, that this FE1+FT 
phase can be considered as a low-temperature bridging structure between the FE1 and FT 
phases since we numerically found that the FEl-to-FEl+FT and FEl+FT-to-FT transi- 
tions are second-order in character, unlike the FT-to-FEl phase transition that displays all 
the expected features of a first-order phase transition (e.g., the toroidal moment suddenly 
disappears at this boundary in favor of a finite value of the polarization, and the FT-to-FEl 
boundary line has a rather large thermal hysteresis: it is typically increased by ~ 100 K 
with respect to the one displayed in Fig. 1 when heating, rather than cooling). 

The last two phases, FE2 and FE3, appearing in Fig. 1 are both ferroelectric and occur 
at intermediate and low temperatures, respectively, for the largest negative values of An. 
In other words, these two states correspond to cases for which the dot, unlike the medium, 
is made of a material that desires to be paraelectric. As a result, the FE2 phase exhibits a 
spontaneous polarization that originates from the alignment along a < 001 > direction of 
dipoles belonging to some specific regions of the medium. More specifically, these regions 
belong to the {001} planes that contain the polarization direction and that do not possess 
any site belonging to the dot (see, e.g, the corresponding inset of Fig. 1 for FE2 showing the 
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homogeneity of dipoles in the top and bottom planes, while the parts of the medium located 
at the right and left sides of the dot do not display any homogeneity for their dipoles). The 
reason behind such unique arrangement, in which the largest dipoles in the medium inten- 
tionally avoid to point towards the dots, is once again the minimization of the depolarizing 
energy. When decreasing the temperature, the FE2 phase transforms into the FE3 state, 
that is associated with a polarization that has now two non- vanishing components along two 
different < 001 > directions. More precisely, the corresponding inset of Fig. 1 reveals that, 
in the FE3 phase, the top and bottom parts of the medium exhibit dipoles that are similar 
in direction than those in the FE2 phase, while dipoles located in the medium at the right 
and left sides of the dot have dipoles aligned along a perpendicular direction. Such dipole 
arrangement arises, once again, from a minimization of the depolarizing energy. Interest- 
ingly, the FE2 and FE3 phases bear resemblance to some states that were experimentally 
found recently in artificially constructed magnets 19 . Moreover, note that we numerically 
found that another solution (but of slightly higher energy than that of the FE3 phase) is 
possible at low temperature for large negative Ak: it consists of the top and bottom parts 
of the medium having dipoles homogeneously aligned a specific < 001 > directions (as in 
the FE2 phase) while the dipoles of the medium located at the right and left sides of the 
dot form a vortex structure. 

We also numerically checked that all the phases of Fig. 1 still occur when varying the size 
of the dot or the size of the whole supercell. Moreover, out of these six states, the FT and 
FE1+FT phases are the soles structures that refine themselves when allowing several AB"03 
dots to be present in a supercell possessing a AB'03 medium. More precisely, Figs 2 show 
that, for large enough positive Ak, neighboring dots have vortices rotating in an opposite 
fashion. In other words, the FT (respectively, the FE1+FT) phase of Fig.l should in fact 
become the AFT (respectively, FE1+AFT) structure of Fig. 2a (respectively, Fig. 2b) when 
dots are close enough from each other (as in dots' arrays). Such novel antiferrotoric phases, 
unlike the FT and FT+FE1 states, allow the dipoles located between two adjacent dots to 
all point along similar directions, which minimizes the corresponding domain wall energy at 
the cost of the electrostatic interaction between neighboring vortices (we numerically found 
that the direct electrostatic interaction between two different toroidal moments should align 
these moments along the same direction - that is, it should lead to a ferrotoric, rather than 
antiferrotoric, order). Furthermore, Fig. 2b reveals that the centers of the vortices in the 
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FE1+AFT phase form a rather unusual lattice (see, e.g., that the centers of the vortices 
have different y-locations for dots that are adjacent along the x-axis). This unique geometry 
originates from the desire of the whole system to maximize the number of its dipoles (in the 
dots) lying along the polarization direction, considering the underlying antiferrotoric order. 

In summary, varying the polarizabilities of the embedded dots and of the medium opens 
a new way for engineering dipole patterns. The fundamental reasons behind such engineer- 
ing are depolarizing-fields effects and medium-driven interactions between dots. Interest- 
ingly, our findings may result in promising applications using negative refraction materials^, 
nanomemory devices, etc.... They should also be relevant to, and thus lead to revisiting the 
analysis of, samples (i) made of composites^; (ii) in which segregation takes places (e.g., 
Ba(Zr,Ti)C>3 for intermediate compositions 22 ); (iii) in which chemically-ordered nanoregions 
form inside a disordered matrix (e.g., relaxors such as Pb(Mg,Nb)Oar^). One can also easily 
imagine that further considering multi-component (e.g., three- component) systems can lead 
to even richer phase diagrams. Furthermore, designing some specific lattice of ferroelectric 
dots (inside some medium) may lead to other exotic states - as similar to the recent finding 
of Ref.— for nanoscale ferromagnetic islands. 
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MRI Grant 0421099 from NSF. S.P. also thanks RFBR grants 01-02-629&05-02-90568_HHC. 



1 J. F. Scott, Ferroelectric Memories (Springer- Verlag, Berlin, 2000). 

2 I. Naumov, L. Bellaiche and H. Fu, Nature 432, 737 (2004). 

3 I. Ponomareva, I. I. Naumov, I. Kornev, Huaxiang Fu, and L. Bellaiche, Phys. Rev. B 72, 
140102(R) (2005); I. Ponomareva, I. I. Naumov, and L. Bellaiche, Phys. Rev. B 72, 214118 
(2005). 

4 H. Fu and L. Bellaiche, Phys. Rev. Lett. 91, 257601 (2003). 

5 V. M. Dubovik, M. A. Martsenyuk and N. M. Martsenyuk, JMMM 145, 211 (1995). 

6 P. S. Neelakantaswamy, B. V. R. Chowdari and A. Rajaratnam, J. Phys. D: Appl. Phys. 16, 
1785 (1983). 



7 



7 Z.-M. Dang, Y. Shen, and C.-W. Nan, Appl. Phys. Lett. 81, 4814 (2002). 

8 G. Lim, J. Osman, and D. R. Tilley, Ferroelectrics 255, 59 (2001). 

9 L. Bellaiche, A. Garcia and D. Vanderbilt, Phys. Rev. Lett. 84, 5427 (2000); Ferroelectrics 266, 
41-56 (2002). 

R. Hemphill, L. Bellaiche, A. Garcia and D. Vanderbilt, Appl. Phys. Lett. 77, 3642 (2000). 

1 L. Bellaiche and D. Vanderbilt, Phys. Rev. B 61, 7877 (2000). 

2 W. Zhong, D. Vanderbilt and K.M. Rabe, Phys. Rev. Lett. 73, 1861 (1994); Phys. Rev. B 52, 
6301 (1995); 

3 C. Bungaro and K.M. Rabe, Phys. Rev. B 65, 224106 (2002). 

4 A. A. Gorbatsevich and Yu. V. Kopaev, Ferroelectrics 161, 321 (1994). 

5 A. Garcia and D. Vanderbilt, in First- Principles Calculations for Ferroelectrics: Fifth Williams- 
burg Workshop, R.E. Cohen, ed. (AIP, Woodbury, New York, 1998), p. 53; Appl. Phys. Lett. 
72, 2981 (1998). 

6 K. Binder, Z. Phys. B 43, 119 (1981). 

7 K. Yu. Guslienko, Appl. Phys. Lett. 78, 3848 (2001). 

8 S. Prosandeev, I. Ponomareva, I. Kornev, I. Naumov, and L. Bellaiche, Phys. Rev. Lett. 96, 
237601 (2006). 

9 L. J. Heyderman, F. Nolting, D. Backes,S. Czekaj, L. Lopez-Diaz, M. Klaui and U. Riidiger, C. 
A. F. Vaz, J. A. C. Bland, R. J. Matelon, U. G. Volkmann, and P. Fischer, Phys. Rev. B 73, 
214429 (2006). 

;0 I. B. Vendik, O. G. Vendik, M. S. Gashinova, Tech. Phys. Lett. 32, 429 (2006). 
:1 C. Huber, M. Treguer-Delapierre, C. Elissalde, F. Weill, and M. Maglione, J. Mater. Chem., 
13, 650 (2003). 

;2 R. Farhi, M. El Marssi, A. Simon, and J. Ravez, Eur. Phys. J. B 9, 599 (1998). 
:3 M. A. Akbas, and P. K. Davies, J. Phys. Chem. Solids 61, 159 (2000). 

' A R. F. Wang, C. Nisoli, R. S. Freitas, J. Li, W. McConville, B. J. Cooley, M. S. Lund, N. Samarth, 
C. Leighton, V. H. Crespi, and P. SchifTer, Nature 439, 303 (2005). 



8 



CAPTIONS 



Fig. 1. Temperature-versus-Afc phase diagram of a 12 x 12 x 12 AB"0 3 dot embedded 
in a AB'0 3 medium within a 16 x 16 x 16 periodic supercell.The positive Ak part of this 
diagram corresponds to a soft ferroelectric dot immersed in a medium that is ferroelectrically 
harder than the dot and that has a decreasing ferroelectric instability as Ak increases. The 
negative Ak part of this diagram corresponds to a dot (having a ferroelectric instability that 
is weaker than those of the medium and that decreases, and then vanishes, as Ak increases in 
magnitude) embedded in a ferroelectrically-soft medium. The lines with symbols represent 
the phases' boundaries. The insets show a (001) cross-section of the dipole configuration in 
the different phases. Specifically, these insets correspond to atomistic calculations with the 
following (Ak, temperature) combination: (-0.0212 a.u., IK), (-0.0212 a.u., 500K), (0.0062 
a.u., IK), (0.0087 a.u., 1 K) and (0.0112 a.u., IK) for the FE3, FE2, FE1, FE1+FT and FT 
phases, respectively. The dot surfaces are indicated via thick continuous lines in these insets. 
The x- and y-axes are chosen along the pseudo-cubic [100] and [010] directions, respectively. 

Fig. 2. (001) cross-sections of the dipole configuration in the AFT and FE1+AFT phases 
in panels a and b, respectively, for four 12 x 12 x 12 AB"03 dots embedded in a AB'03 
medium within a 32 x 32 x 32 periodic supercell. Such cross-sections specifically correspond 
to atomistic calculations for which the (Ak, temperature) combination is (0.0112 a.u., IK) 
and (0.0087 a.u., IK) for the AFT and FE1+AFT states, respectively. The dots surfaces are 
represented by thick continuous lines. The x- and y-axes are chosen along the pseudo-cubic 
[100] and [010] directions, respectively. 
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